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ABSTRACT 

We examined the Wilson-Bappu effect, a relationship between the absolute mag- 
nitude of the star, My, and the logarithm of the Ca n emission width, Wo, over the 
largest My range to date, +13 to -5, covering M-dwarfs to type la supergiants. We used 
an extensive literature, the latest Hipparcos reduction, data from two globular clusters, 
and new observations from Apache Point Observatory to compile a sample that allowed 
us to study the effect of [Fe/H] on the Wilson-Bappu relationship. Our results include 
reporting the deviations from linearity and demonstrating that the Wilson-Bappu rela- 
tionship is insensitive to metallicity. 

1. Introduction 



O.C. IWilson fe Vainu Bappul (|19571 ) published a remarkable linear correlation between the 
logarithm of the Ca II K emission full width at half maximum (log Wo) and the absolute visual 
magnitude of the star, My. They presented that this holds for non- variable cool stars below type 
F, regardless of spectral type, [Fe/H], or intensity of the emission line. This correlation is useful 
to derive an My for cool sta rs too distant for accurate paralla xes. This relationship can also be 
found in other spectral lines, lAmbrusterl (jC. W. and Pettersenl ) were the first to publish a similar 
relationship between My and the with of the Mg n k line. 

The Wilson-Bappu (hereafter WB) effect has ai ded i n studies of var i able stars like T Tauri 
stars as well as main sequence stars (see iKuhi I ()1965l ) and lBaliunas et al.l (|2005l )). The WB effect 
ha s been invaluab l e for ou r abili ty to under s tand c hromosph eres and he ating, see seminal work 
by iReimersI ()1973l ) , iNeckell ()1974l ) and lAyresI (| 19791 ) . IReimersI and iNeckell derived a way to relate 



Wq to stellar parameters: W ~ g T [Fe/H] . lAyresI was able to establish that the WB 



effect is largely a result of the surface gravity of the star. The Ca n line is so optically thick that 
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the increase in width corresponds to an increase in column density in the chromosphere, 
relationship between column density and Wq is predicted to be insensitive to [Fe/H]. 



This 



Accurate Hipparcos parallaxes were used for the first time I Waller stein et al.l (|l999l ) to study 
the WB effect for a larg e number of sta rs with accurate My. T heir log Wq were taken from the 
large list published by IWilsonl (|1976l ) . IWallerstein et al.1 (|1999i ) were the first to use Hipparcos 
parallaxes and noticed that their small sample of supergiant stars showed unusually large log Wq 
and deviated from a linear trend. 



Wilson &: Vainu Bappu (jl957l ) foudn that the WB effect has no [Fe/H] dependence. How- 
ever, they had relatively few stars compared to recent studies of the WB effect and there are 
few metal-poor stars i n the solar neighborhood. A chromospheric study of metal-poor giants, 
(jDupree &: Smith! Il995l ) .hereafter DS95, demonstrated that extremely metal-poor stars ([ Fe/H] < 
-2.0) h ave a diminished Wq and did not follow the WB relationship previously established by I Wilson 
(j 19761 ). The field giants they studied were too far away to have accurate parallaxes, therefore they 
estimated their lumino sity by fitting e ach st ar's colors to isochrones of M92, an extremely metal- 
poor globular cluster. ICacciari et al.l (120041) compar e d sta rs they studied in the globular cluster 
NGC2808 to the study of field stars by iPace et al.l (|2003l ) and came to a similar conclusion as 
Dupree fc Smith! (119951 ). where more metal poor stars lie above t he linear fit deriv ed from more 
metal rich stars. However, their result was much more subtle then lDupree fc Smith 



We e xpand the work by IWallerstein et al.l (|1999l ) and DS95 by using the latest Hipparcos 
reduction (jvan Leeuwenl 120071 ) and constraining our sample to stars with parallaxes greater than 
3 times the probable error (Tl/an > 3)0 instead of (H/au > 5). We have included all available 
measurements of Wq from the literature, additional metal-poor stars observed at the Apache Point 
Observatory, and Ca n widths of stars from the globular clusters NGC2808 and oj Cen. 



The Database 



Parallax is the most reliable and direct method of deriving My. Hippa rcos constitute the 
largest sample and most reliable parallaxes, so we use d Hipparcos p arallaxes (jvan Leeuwenl 120071 ) 
for stars that meet our criteria of H/cru > 3. The Ivan Leeuwenl reduction of Hippa rcos data 
reduced the uncertainty in the parallaxes by a factor of 2 from iPerrvman fc ESAl (|1997l ) , allowing 
us to study many more stars with acceptable parallaxes. 

To investigate the WB effect and its dependence on [Fe/H], we have assembled data from a 
list of previously published Ca II emission widths, Wq. Wq is defined as the width of the line, 
W, having the instrumental profile removed. See figure [2] for an example of how to measure W. 



1 When studying individual stars, using parallaxes that are less than 5 times their uncertainty can introduce large 
errors and unreliable results. For a substantial set of data from which the general behavior of the sample is to be 
understood, a less strict limit provides sufficient accuracy for the purpose at hand. 
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Wo values came from the catalog s by iGlebocki et al.l (119801 ). lLastennet Freire Ferrerol (|1994l ) 
as well as IZarro Rodeersl (jl983T ). IWallerstein et al.l (jl999h . and IPace et al.l (|2003h . For the M- 
dwarfs in our samp l e, we used Ca n K full width at half maximum (FWHM) measurements from 



Rauscher &: Marcvl (|2006l ). The sample of Wq values still lacks a sufficient number of metal-poor 



stars, so we included the data by DS95 (shown in Table [T]). In the case of a star with multiple 
observations, the most recent value of Wq was used. For consistency, we used only values of 
Wq. Most of these publications used Wq, but some published values of W. The instrumental 
profiles of the spectrographs used were quoted in the paper and we subtracted the value from their 
published W to obtain Wo . We compared the different references and methods to look for systematic 
differences in Wo that would affect our study. There was no dif ference in the Wilson-Bappu e ffect 
between the different references. We chose to use the catalog bv lCavrel de Strobel et al.l (|200ll ) for 
our values of [Fe/H] for consistency. 

To stud y other main sequence star s, we analyzed the spectra from the McDonald Observa- 
tory archive (jAllende Prieto et al.l 12004 ) and measured Wo ourselves. The spectra from McDon- 
ald were already reduce d and flux calibrated, we only measured the width of the Ca II K line; 



sec 



Allende Prieto et al.l (J2004J) for the full list of stars. We also expanded the sample of metal 



poor stars with observatons with the ARC echelle spectrograph (ARCES) on the 3.5m telescope 
at Apache Point Observatory (APO). Observations were made with the standard 1.6"x 3.2"slit. 
Our spectra have a resolution of 31,500, giving us a spectral profile of 7km/s. ARCES covers a 
wide spectral range of 3,500A-10,000A. We also observed some metal-poor giants at Kitt Peak 
National Observatory on Oct 14-17, 2000. The resolution of the spectrograph used is 18,000. 
We reduced our spectra with IRAF using a data reduction guide, provided on the APO home- 
page (www.apo.nmsu.edu/arc35m/Instruments/ARCES/images/echelle_data_reduction_guide.pdf) 
and the reduction script writen by John Barentine (http:/ /www. apo.nmsu.edu/arc35m/Instruments/ARCES/echel 
We took extreme care to get the wavelength calibration correct, since it is the AA that determines 
W . Our measurements, observation dates, and exposure times are presented in Table [2 

DS95 is of particular interest due to their result of an [Fe/H] dependent WB effect. However, 
we choose to only include Wo measurements that had a parallax above our 3a criteria. Nine of the 
twenty- four stars with Wo have accurate parallaxes in Hipparcos. However, as we show in Table 
[T]and figure [H these stars' My are on average brighter than My derived from the 2007 parallaxes. 
DS95 used M92 to estimate the My of all their stars, which have a variety of [Fe/H]. M92 has a 
luminous and relatively blue RGB due to its low metallicity. When fitting stars of lower metallicity 
to a bright blue RGB, the stars will have a higher estimated My than the star's true My. We 
decided to leave out the stars that do not have parallaxes greater than 3a due to the inconsistency 
of using the RGB to estimate My. Many of these metal poor stars are lost to us since we cannot 
trust the method used to estimate My. 



Pace et al. (j2003l ) and Rauscher Sz Marcy ( 2006 ) used Gaussian fitting programs to determine 



W because the stars they observed have strong emission profiles. However, our APO sample is very 
metal-poor and have small Ca n K amplitudes is not easily fit with a gaussian profile. Therefore, we 
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chose to measure the profile emission width using the method of I Wilson! (119761). as shown in Figure 



H T his method was also used on the full sample of spectra from McDonald (jAllende Prieto et al 



2004! ) . Allende Prieto used the Harlan J. Smith 2.7m telescope and the ESO 1.52m telescope to 
survey nearby stars with a spectral resolution of ~ 50,000. We compared the WB effect for stars 
measured by gaussian fitting vs those stars measured using the FWHM and found that there are 
no systematic differences in Wq space. For the spectra we obtained, we measured the same value of 
W (within a tolerance) using the two methods. The tolerance becomes negligble after taking the 
logarithm. 

We include studies of globular clusters NGC 2808 and u Cen to expand our sa mple of metal 



poor s tars, allowing us to compare the WB effect among multiple stellar populations. ICacciari et al 



(|200J) provided measurements of log W n and V and (m-M v ) of RBG stars in NGC 2808. [Fe/H] 
values of these stars were published bv lCarretta et al.1 (120041). The measurements are shown in Table 
El Omega Cen has distance modulus of -14.07 ( Harris! 1996 . the 2010 edition). The V and [Fe/H 



for st ars ROA159, 256, 238, 523, WFI321293, and WFI 140419 w e re pu blished by IVievtes et al 



(2011 



The identifier ROA stars are from the catalog by IWoolleyi (|1966l ) and the stars identified 
by W FI were observed wit h the Wide Field Imaging camera on the 2.2m ESO-MPI telescope in 
Chile (jPancino et al.l l2000i ) . We obtained the spectra used inthese studies and measured Wo of 
these stars using the same method as APO and McDonald, described above. The uj Cen stars are 
presented in Tabled! 



3. Fitting the Wilson-Bappu effect 



The WB effect was thought to be linear in log Wq until lWallerstein et al.l (119991 ) demonstrated 
that the brightest stars hav e smaller luminosities than their Ca H widths would suggest. M-dwarfs 
showed nonlinearity in Mg il lElgaroy et al. I (|l990l ) and Rausc her Marcvl (120061 ) observed unusual 
nonlinearity in Ca n FWHM. We used the stars in Rauscher &z Marcy! ( 2006 ) to compare against 
the other stars in our sample. Rauscher Marcv ( 20061 ) noted that the Ca n FWHM of those stars 
were within a few km/s of the intrinsic instumental profile. 

Our full sample shown in Figure [3] with a weighted least squares fit of the linear part of the 
plot (luminosity classes II, III, IV, and bright V's). Eqn ([T]) is the line of best fit for stars that 
fall between 1.4 < log Wo < 2.0 where the errors quoted are the standard deviation of the fit. 
It is evident that the highest and low est luminosity class star s do not follow the linear trend of 
the giants as previously pointed out by IWallerstein et al.l (|1999l ). As a group, the supergiants (and 
perhaps some misclassified luminosity class II stars), follow a different trend than the rest of the 
giants. There seems to be a step discontinuity at the transition between the M-dwarfs and the 
K-dwarfs (log Wo ~ 1.4) where there is a spread in My- Accounting for discrepancies in the 
dwarfs and supergiants is beyond the scope of this study and no higher order polynomial fit was 
attempted. 
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My = -14.7(±0.29) log W + 27.0(±0.52) 



(1) 



3.1. Low Mass Main Sequence Stars 

Main Sequence stars have a strong and well established relationship between surface gravity, 
luminosity, and effective temperature. Therefore, we wanted to comapre the expected and observed 
low luminosit y WB effect. Using the log g, T e jf, and My from the Zero- Age Main- Sequence 



i sochr ones b y iGirardi et al.l (|2000l ) and the relationship (Wq ~ g T ) established by iReimers 



(j 19731 ) and iNeckell (|1974l ) we can compare expected and observed WB effect. The expe cted and 
observ ed low luminosity WB effect is shown in Figure [H The WB effect expected from IReimers 
(j 19731 ) should be linear down to My =13. 



Using the conclusions of lAvresI (|1979l ). we can speculate that the column depth of M-dwarf 
chromospheres is a lot smaller, leading to a thinner Wq that we observe. We can use the same 
argument to say that supergiants have a much thicker chromosphere then models predict, leading 
to a larger Wq than predicted. 



3.2. Metallicity Dependence 

Because of the large sample and well-behaved linear trend, the luminosity class II, III, IV, and 
V stars are best suited to examining the effect [Fe/H] has. Figure [5] shows residuals from Eqn 
([I]), A (log Wq), as a function of [Fe/H]. The averages of A log Wo for the metal-poor stars and 
the metal-rich stars are —0.034 ± 0.017 and —0.048 ± 0.042, respectively. We conclude that on 
average, the stars with [Fe/H] < -1.0 do not deviate significantly from the linear fit. Comparing 
A (log Wo) between the metal-poor stars and metal-rich stars shows that they are drawn from the 
same distribution. A Kolmogorov-Smirnov test between the metal-poor and metal-rich samples 
yields a value of 0.27 with a significance of 0.0005, confirming the null hypothesis that the two 
groups are drawn from the same parent distribution. 



A log W = 0.027(±0.004) + 0.0005(±0.004)[.Fe/.ff] 



(2) 



The lack of [Fe/H] trend disagrees with the conclusions of DS95. The latest Hipparcos My 
differ from the values that they derived using M92 (see Figure [1]). Using only My from Hipparcos 
parallaxes, the Pearson's r correlation c oefficient gives a valu e of 0.09, indicating that A log Wo is 
uncorrelated with [Fe/H ] , The study bv lCacciari et al.l (|2004l ) compare stars from NGC 2808 with 
the study by lPace et al.1 (120031 s ) noting tha t stars in NGC 2808 had a wider spread and tended to be 
brighter than the fit by lPace et al.1 (|2003l ) would suggest. However, this could just be a sampling 
effect. Our study shows conclusively that any dependence on [Fe/H] is has a much smaller effect 
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on the WB effect than log g. 

We preserve equation [1] as it stands, adding no term for [Fe/H]. More accurate My and a larger 
sample of log Wq will be required for a better WB relation. A re-analysis of the [Fe/H] dependence 
will be possible after the Gaia satellite measures parallax better than 0.3 mas out to V = 20 and 
better than Wfias at V=15 magnitude. By comparison, Hipparcos measured parallaxes down to 
5mas at a limiting magnitude around V=10. 



4. Conclusion 

We present the largest sample of stars to demonstrate the Wilson-Bappu effect over the largest 
range of My to date. The main results of this paper are shown in Figs [3] and Our large sample 
shows that there is a limited range where the WB effect is linear between +8 < My < —3 and 
1.35 < log Wq < 2.0. Stars deviate from Eq 1 above and below these ranges. To our knowledge, 
there is no theory of the Ca n emission width that predicts that the relationship should be linear 
from 13 < My < —6 or that a higher-order polynomial would be preferable, the issue of what 
curve should be drawn through the points is not important. What is important is what we can 
learn from devitation from the luminosity classes. Ayres (1979) was one of the last developments 
on width-luminosity relationships like the Wilson-Bappu effect making arguments that increasing 
column density leads to an increase in the Ca II width. We can come to the conclusion that the 
column density of M-dwarf chromospheres are lower than expected while the supergiants have a 
larger chromospheric column density than expected. 

Modern instrumentation such as the Hipparcos satellite and CCD detectors on ground-based 
telescopes have permitted the WB effect to be extended to about My=+13. The value of Figure 
3. is the recognition of stars that deviate significantly from the mean so that so that their chro- 
mospheres may be investigated further in detail. Variations in the intensity of the Ca n lines in 
G and K dwarfs have yielded a great deal of informa tion on their rotation and solar-type cycles of 



varying chromospheric activity iBaliunas et al.l (|2005l ) . Our data provides an observing list of stars 



whose monitoring should produce very interesting results, especially for giants and supergiants of 
types G, K, and M on their rotations and chromospheric activity variations. 

The improved 2007 Hipparcos reduction has allowed us to determine My of stars farther away, 
including more metal-poor stars. High resolution spectra of metal-poor globular cluster giants 
allowed us to explore the dependence of [Fe/H] to answer some fundamental questions about the 
WB effect: Is metallicity the cause of the 0.3 magnitude spread, or is it responsible for the deviations 
of the high and low luminosit y trends?. Figure 4 demonstra tes clearly the lack of depende nce on 



[Fe/H ], as initially asserted by I Wilson k, Vainu Bappul (Il957l ) and theoretically developed by I Ayres 



fliaza ). 
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Table 1. Metal- Poor Giants from 



Dupree fc Smith 



1 fi tin "i - 1 ti Q7 1 
XUCllLlll"! 


log 


V'O 


MQ9 Mir 


-1-1/ on 


H 1 T^T~\ OTPriG A/Ytt- 
111U Ucll LUo 1VJ. y 


fFp /Hi 


HD6268 


1 


78 


-1.2 






-2.36 


HD6833 


1 


83 


-0.9 


8.5 


0.31± 0.25 


-1.04 


HD8724 


1 


92 


-1.3 


5.5 


0.47±0.43 


-1.83 


HD25532 


1 


90 


-04 


4.1 


1.36±0.53 


-1.26 


HD26297 


1 


83 


-0.9 






-1.68 


HD29574 


1 


87 


-2.6 






-1.88 


HD63791 


1 


78 


-0.6 






-1.81 


HD103036 


2 


00 


-2.2 






-1.78 


HD110184 


1 


77 


-2.3 






-2.44 


HD110281 


1 


99 


-2.0 






-1.56 


HD118055 


1 


87 


-1.9 






-1.76 


HD122563 


1 


79 


-1.5 


12.1 


-0.55±0.18 


-2.74 


HD122956 


1 


80 


-0.9 


5.4 


-0.12±0.41 


1.63 


HIP70199 


1 


90 


-2.3 






-1.99 


HIP73960 


1 


86 


-1.7 






-1.37 


HD165195 


1 


84 


-1.8 






-2.24 


HD166161 


1 


85 


-1.2 


5.4 


1.56±0.40 


-2.24 


HD175305 


1 


74 


1.2 


17.8 


1.35±0.12 


-1.16 


HD 184266 


1 


96 


-1.5 


7.0 


1.19±0.31 


-1.44 


HD187111 


1 


81 


-1.9 






-1.54 


HD204543 


1 


92 


-0.3 






-1.84 


HD216143 


1 


85 


-1.5 






-2.25 


HD221170 


1 


81 


-1.96 


4.3 


0.18±0.51 


-2.19 


HD232078 


1 


90 


-2.0 






-1.54 



-10- 



Table 2. APO Observations 



Identifier 


My 


log 


Wn 


[Fe/Hl 


Date Obs(UT) 


Exposure time 


HD6755 


2.14±0.25 


1 


66 


-1.58 


2010-11-21 1:33:36 


lOmin x2 












2010-08-21 11:07:27 


5 min 


HD6833 


0.31±0.25 


1 


73 


-1.04 


2010-08-21 10:59:38 


3 min 


HD25532 


1.36±0.52 


1 


79 


-1.26 


2010-04-12 


5min 


HD37828f. 


-0.98±0.40 


1 


84 


-1.42 


2000-10-14 


2hrs 


HD73394 


-0.34±0.69 


1 


80 


-1.40 


2011-04-12 03:46:16 


20 min 
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1 


84 
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1 


85 
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1 


89 


-1.85 


2000-10-15 


lhr 
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1 


66 


-2.19 


2010-08-21 10:09:10 


5min 



Table 3. \ogW and [Fe/H] for NGC2808, Identifiers from (|Cacciari et all 12004 ) 



ID No. 


log W 


My 


[Fe/H] 


37872 


1.86 


-1.940 


-1.06 


46099 


1.93 


-1.849 


-1.19 


46580 


1.89 


-1.900 


-1.03 


47606 


1.93 


-2.154 


-1.14 


48889 


1.91 


-2.249 


-1.22 


50119 


1.87 


-1.704 


-1.21 


50761 


2.06 


-2.200 


-0.81 


51454 


1.90 


-2.144 


-1.11 


51499 


1.99 


-2.155 


-1.24 
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Fig. 1. — The graphical representation of Table [TJ demonstrating that the My estimated by DS95 
do not accuately derive the same My from Hipparcos parallaxes. The DS95 absolute magnitudes 
are systematically brighter than My by Hipparcos. The solid line shows when both values of My 
agree. This only includes stars with parallax greater than 3 times the probably error. 
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Fig. 2. — The Ca n K emission profile of HD64394 observed from McDonald. W is measured as 
the wavelength difference of the half maximum of the blue and red side of the emission profile at 
the bottom of the K line. W is measured in km/s. The instrument profile is then subtracted to 
get Wq, modern instruments have an instrumental profile less than 10km/ s. 
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Fig. 3. — The Wilson-Bappu effect shown is broken down by luminosity class. The dotted 
fit is defined in Equation 1 and is the weighted least squares fit of the linear region between 
1.4 < log Wq < 2.0. The luminosity class I stars seem to follow a trend that takes them below 
what Eqn (1) would predict. The M-dwarfs also seem to fall above the fit, having smaller Wq than 
predicted. 
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Fig. 4. — The expected Wilson-Bappu effect is shown by the black lines, parameters derived from 
theoretical isochrones. The data shows that the M-dwarfs have significantly diminished Ca II widths 
than predicted. This would indicate a smaller column density in the chromosphere than the stellar 
models indicate. 
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Fig. 5. — The deviations of stars from Eqn (JTJ) (log Wq) are shown as a function of [Fe/H]. The 
weighted least squares fit is given in equation (2). This plot shows that [Fe/H] has no effect on the 
deviations from linearity of the giants. This confirms the [Fe/H] dependence predicted by Reimers, 
Neckel, and Ayres. 
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Table 4. The u Cen Sample 



Identifier 


log W 


V magnitude 


[Fe/H] 


ROA500 


1.91 


13.12 


-0.53 


ROA517 


1.92 


13.10 


-0.54 


WFI139267 


1.73 


13.09 


-0.79 


WFI263340 


1.90 


13.62 


-0.63 


WFI305654 


1.87 


13.38 


-0.65 


ROA159 


1.89 


11.98 


-1.72 


ROA256 


1.95 


12.28 


-1.71 


ROA238 


1.84 


12.44 


-1.80 


ROA523 


1.82 


13.39 


-0.65 


WFI321293 


1.88 


13.69 


-0.72 


WFI140419 


1.86 


13.49 


-0.68 



